Ion and neutral molecules in the W43-MM1(G30.79 FIR 10) infalling clump 

Paulo C. Cortes^ 

National Radio Astronomy Observatory - Joint ALMA Office, Alonso de Cordova 3107, 

Vitacura, Santiago, Chile 

pcortesOalma . cl 
Received ; accepted 



-2- 
ABSTRACT 

The high mass star forming clump W43-MM1 has been mapped in N2H*(4 — > 3), 
C^^O(3 2), SiO(8 7), and in a single pointing in DCO+(5 4) towards the 
center of the clump. Column densities from these observations as well as previous 
HCO+(4 ^ 3), Hi3CO+(4 ^ 3), HCN(4 ^ 3), H"CN(4 ^ 3), and CS(7 ^ 6) data, 
have been derived using the RADEX code, results later used to derived chemical abun- 
dances at selected points in the MMl main axis. Comparing with chemical models, 
we estimate an evolutionary age of 10'* years for a remarkable warm hot core inside 
MMl. We also proposed that the dust temperature derived from SED fitting in MMl 
is not representative of the gas temperature deep inside the clump as dust emission 
may have become optically thick. By deriving a deuterium fractionation of 1.2 x 10^, 
we estimate an electron fraction of X(e) = 6.5 x 10"^. Thus, the coupling between the 
neutral gas and the magnetic field is estimated by computing the ambipolar diffusion 
Reynolds number i?;;, = 18 and the wave coupling number W = 110. Considering that 
the inf ailing speed is slightly supersonic (M=l.l) but sub-alfvenic, we conclude that 
the MMl clump has recently or is in the process of decoupling the field from the neu- 
tral fluid. Thus, the MMl clump appears to be in an intermediate stage of evolution 
in which a hot core has developed while the envelope is still infalling and not fully 
decoupled from the ambient magnetic field. 

Subject headings: ISM: Abundances, ISM: Magnetic Fields, ISM: clouds, ISM: Kinematics 
and dynamics 
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1. Introduction 



High mass star formation involves the study of gas masses many orders of magnitude larger 
than low mass stars increasing the complexity as the process is highly energetic and dynamic. 
The distances at which the regions of interest are located are usually kilo-parsecs away from the 
sun, complicating the use of single dish telescopes as current dish sizes do not produce enough 
angular resolutions to resolve individual cores. Thus, the emission is often diluted within the 
beam complicating the interpretation as the high mass cores inside clumps are usually shrouded 
by optically thick envelopes, even with dust emission. We know that massive stars (>10Mq) 
form in giant molecular clouds several orders of magnitude larger than the low-mass star forming 
clouds. As it is still not clear how these giant ensembles of gas came to be, we know that their 
location are correlated with the spiral arms of our galaxy. These massive complexes are known 
to harbor cores with usually several sources as massive stars are born in groups or associations. 
Once a massive star turns on it \ yill quickly perturb and ion i ze the surrounding envelope creating 



what it is known as a hot core (|Wood & Churchwell 



19891 : 



HoareetaL 



20070. The intense 



radiation coming from these young sources, will evaporate C and O bearing molecules from 



ice ma ntles surrounding dust grains, changing dramatically the chemical diversity around the 



star(s) (iHerbst & van Dishoeck 



2009'). Even though, the inventory of complex molecules in hot 



cores have been growing consistently in the past 10 years, still chemical and dynamical models 
have problems tr ying to reproduce a bundances and explain the chemical differentiation seen in 



high mass cores (IGarrod et al. 



20081). Thus, probing the densest regions becomes paramount to 



understand their physical and chemical conditions. On the otherhand, the complete set of physical 
parameters have to be considered if we want to fully understand the process. Magnetic fields are 
known to be ubiquitous in the ISM and giant molecular clouds are no exception. Strong magnetic 
field strengths have been measured fro m CN Zeeman splitting towa rds high mass star forming 



regions to be in the order of milligauss (|Crutcher et al. 



1996 . 



19991). Also, observations done 



from polarized dust emission find ordered field structure with strength estimations around the 
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same values (ILailllOOlUCortes et alJl2005l : ICortes & Crutchedl2006l : ICortes et al.ll2008l) . Thus, 
incorporating the magnetic field along with the set of parameters used to study the star formation 
process is critical. 

In this paper we present new observations of DC0"^(5 4), N2H'^(4 3), and SiO(8 — > 7) 



as well as re-processecLl HC0+(4 ^ 3), H" CQ+(4 ^ 3) 



CS(7 6) observations done previously by 



Cortes et al. 



HCN (4 ^ 3), H"CN(4 ^ 3), and 



(120101) and hereafter PCI, from the 



high mass star forming clump W43-MM1 (also known as G30.79 FIR 10). The W43-M M1 is a 



Liszt 



arge high mass star forming clump located within the W43 region near / = 31°, b = 0°. 
(19951) observed the W43 molecular complex in HCO^ and '^CO finding a series of several rings 



and shells in the den se molecular gas, which they attributed to be a product of star formation. 



Moonev et al. 



(119951 ) made the first dust continuum observations of this source at 1.3 mm using 



the IRAM 30-meter telescope. They detected a total flux of 13.6 Jy positioning MMl at an 



interface with the extended H II region in the W43 complex. Several H2O masers (|Cesaroni et al 



19881) have also been detected within an arcsecond from the dust peak found in the 1.3 mm data. 
No centimeter radio-continuum emission seems to be associated with MMl, which suggests that 
the source is in an early stage of evolution. Additional continuum mapping has been done by 



Motte et al. 



(|2003l) at 1.3 mm and 350 fim with the IRAM 30-m and CSO telescope^l respectively. 
They determined a v/.v,. = 98.8 km s ' from H^^CO'^ and a dust temp erature of T^^^g^- ~ 19 K, b y 



Cortes & Crutchei (200 



3) 



fitting a grey body model to the spectral energy density distribution, 
made interferometric observations of polarized dust emission with BIMA, finding an ordered 
pattern for the field and estimating and deriving a total mass of about 3000 Mofor the MMl clump. 



'By re-processing, we mean using the data cubes obtained by PCI to produce maps and calcu- 
late column densities with RADEX (see section [2] which have not been done previously 

^IRAM stands for Institut de Radioastronomie Millimetrique, and CSO stands for Caltech Sub- 
millimeter Observatory 
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Finally, PCI studied the gravitational equilibrium of this source finding compelling evidence for 
infalling m otions and refining a prev ious estimation of the magnetic field strength in the plane 



of the sky (|Cortes & Crutchei 



2006|) to B = 855 fiG. PCI also produced several maps from the 



previously mentioned molecular species. From this set of observations, we have derived chemical 
abundances as well as spatial distributions of molecules, which we use to estimate the electron 
fraction and the coupling parameters between the magnetic and the neutral gas. The paper is 
organized as follows. Sect. 1 is the introduction. Sect. 2 presents the observation configuration. 
Sect. 3 the observational results, while in Sect. 4 we discuss the chemical abundances, the spatial 
differentiation, the fractional ionization, and the ion-neutral coupling. Finally, Sect. 5 presents the 
summary and conclusions. 



2. Observation Configuration 

2.1. ASTE Observations 

The W43-MM1 clump was observed during July 2010 using the Atacama Sub-millimeter 
elescope Ex periment (ASTE) from the National Astronomical Observatory of Japan (NAOJ) 



(IKohno 



2005|) . The telescope is located at Pampa la bola in the Chilean Andes plateau reserve for 
Astronomical research at 4900 meters of altitude. ASTE is a 10 m diameter antenna equipped with 
a 345 GHz double side band SIS-mixer receiver. We observed DC0'^(5 4) and N2H+(4 3) 
by tuning the CAT345 receiver to the appropriate frequencies (see Table [T] for frequencies and 
map sizes) giving a beam size of ~ 20", and a velocity resolution of 0.1 km s"^. The MAC (a 
XF-type digital spectro-correlator) was set to high resolution mode with a total bandwidth of 
128 MHz. The pointing accuracy was in the order of 2", with VY-CMa used as the pointing 
source. The observations were done by performing single pointing (position switching) with 
an off position carefully selected to avoid contamination to the on position. Also, a small 
60" X 60" map for N2H^ was obtained centered at the reference position. The observations 
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were done under excellent weather conditions (precipitable water vapor or <PWV> 0.5 mm 
and wind speed under 3 m/s). Our reference position was {a, 6) = (18M7' "46.9% -l°54'29.r0 



(J2000), wh i ch coi n cides with the peak dust emission repor t ed and used by 



Moonev et al 



Motte et al. 



(l2003h : 



Cortes & Crutched (120061) : 



Cortes et al. 



(Il995h : 



(|2010l) : while the OFF position used 



is {a, 6) = (18''49'^46.7^ -2°00'52.0") Initial data reduction and calibration was done using 
the NEWSTAR package, and the calibrated data were later exported into our own software for 
analysis and plotting. 



2.2. APEX observations 

Observations were performed during the first week of August 2008 using the Swedish 
Heterody ne Facility Instrum ent (SHFI) mounted on the Atacama Pathfinder Experiment telescope 



(APEX) dOusten et al 



120061) ■ located at llano de Chajnantor also in the Chilean Andes plateau 
reserve for Astronomical research at. The SHFI was tuned to 347.33 GHz to detect the 
SiO(7 = 8 — > 7) molecular transition. The spectrometer was set up to 8192 channels with a 
resolution of 0.1 kms ^ beam efliciency is 77=0.73+0.07 as measured by the APEX staff", with 
a pointing accuracy better than 2" and a beam size of 19". The APEX wobbler was not used 
for this observations. Jupiter and R-Aql were used as intensity and pointing calibrators where 
the observations were calibrated by the usual chopper-wheel method. The observations were 
done in raster mode with spacings of 15" from the same reference position used for the ASTE 
observations. The initial data reduction was done with the GILDAS-CLASS reduction package 
and the final analysis with our own software tools. 



Table 1 . Observing parameters for molecular lines observations towards W43-MM 1 . 



Line 


Transition 


Frequency 


Beam Size 


Area mapped 


Peak 


Line- width 


/Tdv 


cr 


Telescope 






[GHz] 


[arcsec] 




[K] 


[kms-i] 


[Kkm/s] 


[K] 




DCO+ 


(/ = 5 ^ 4) 


360.1698810 


20.9 


single pointing 


0.044 


2.9 


0.11 


0.01 


ASTE 


N2H+ 


(y = 4 ^ 3) 


372.6725090 


20.3 


60" X 60" 


2.1 


3.7 


7.3 


0.08 


ASTE 


SiO 


(7 = 8^7) 


347.3305786 


21.7 


60" X 60" 


0.6 


n/a 


n/a 


0.24 


APEX 



I 

-J 
I 
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3. Results 



3.1. Neutrals Observations 



In PCI, we reported a single 30" x 30" averaged spectrum, centered at the MMl reference position, 
for C^^0(3 2). Even though we mapped a larger area of 150" x 120" sampled every 15", we 
did not present the complete set of observations in that work. Here, we present the integrated 
intensity map in Figure [TJ The map follows the elongated structure seen from the '^C0(3 2) 
map presented in PCI (see Figure 2 there). Two clear peaks are seen in the C'^0 map, one is 
coincident with the reference position and the other is located at (a, 5)=(-120", -45"). The first 
peak agrees with the peak in the dust emission, HCN(4 3), HC0+(4 3), and CS(7 6) 
presented in PCI. The second and strongest peak is interesting as it is located at the edge of the 
map following the ^^C0(3 — > 2) peak where it appears to be anti-correlated with the dust emission 
i.e. the dust e mission is minimum where CO is maximum. According to the continuum maps of 



MMl done by 



Motte et al. 



(|2003|) 1.3 mm and 350 fim, the dust emission decreases an order of 
magnitude with respect to the dust-peak at this position. As previously suggested by PCI, this 
area seems to be at an interfa ce with the H II reg i on produced by a large cluster of massive stars 



located at the center of W43 (Wilson et al 



I970I : 



Smith et al. 



19781). Thus, it it possible that we 



are seeing part of a PDR at that area. 

We observed an area of 60" x 60" in SiO(8 — > 7) towards the center of the MMl clump (see 
Figure O. Large amount of data were discarded as the baselines showed evident ripples which 
we could not remove. There is no conclusive evidence of SiO emission from these data, only an 
upper limit of 0.6 K around the systemic velocity of MMl is obtained with an rms noise of 0.24 K 
which is less than 3cr. Thus, we are not able to confirm the marginal evidence find by PCI from 
the HCO^ and CS line- wings. However, new un- published SMa| interferometric observations 



^ SMA stands for Sub Millimeter Array 
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dSridharan to be submitted 



201 1|) found evidence for outflow emission towards MMl as the 
interferometer filters the envelope probing angular scales closer to the center of the core. Thus, it 
is likely that the outflow is young and deeply embedded into the clump with its emission diluted 



within the APEX beam. 



3.2. Ion observations 

Single pointing observations were obtained towards the W43-MM1 core in order to detect the 
DC0^(5 4) molecular transition. Figure [3] shows the spectrum obtained and its corresponding 
Gaussian fit (in red); Table [U summarizes the main molecular line parameters. With a single peak 
of 44 mK over a noise level of cr = 10 mK, the line intensity corresponds to a 4cr significance 
observation, which is clearly a detection. The data were processed by keeping only spectra with 
stable baselines using order 1 polynomial fits. The spectrum was binned every 10 channels, giving 
a velocity resolution of 1 kms ^ The DCO^ molecule has been detecte d in its lower transitio n 



towards a number o f star forming cores with similar results to our own dWootten et al 



1979 



Caselli et al. 



2002a 



b|). The line profile is clearly Gaussian, which suggests that the line is likely 
optically thin. Even though it is expected to find high abundances of DCO^ in early stages of 
massive cores, the abundances are unlikely to be large enough to make the line optically thick. 
By looking at Figure 1 we also noticed a strong line besides the DCO*(5-4) detection (with a 
Gaussian fit in blue). The corresponding line center frequency from the best fit is 360.1847 GHz, 
it is currently un-identified and its nature will be discussed in section |4~n 

We observed N2H+(4 —> 3) rotational transition towards the center of MMl. A 60" x 60" 
integrated emission map is presented in Figure IH The map is centered at the reference position 
and the emission was sampled every 10" in both a and 6. Also, a short 5 minutes scan was done 
at the center to increase the signal to noise in the spectrum. The spectrum, shown in Figure |5] with 
its corresponding Gaussian fit, was binned every 10 channels giving a velocity resolution of 1 
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km s ^ From the fit, we obtained a peak intensity of 2.1 K with a mis noise of cr = 80 mK. No 
hyperfine structure is seen as is likely blended due to a large line- width of Av = 7.5 km s~^seen 



in the emission, wh ich has also been seen by others (iFriesen et al. 



Fontani et al. 



2010L for N2H+(4 ^ 3) and 



20061 for N2H^(3 — > 2)). The N2H^ center spectrum seems to be consistent with 
an optically thick line centered at v =97.14 km s \ which is blue-shifted from the core v/vr and 
in agreement wit h infalling evidenc e presented in PCI and by others for similar high mass star 



forming regions (|Fuller et al. 



|2005|) . The N2H^ integrated intensity map appears as a filament 
from South-East to North- West, which is consistent with the ^^C0(3 2), HCN(4 3), and 
HC0^(4 — > 3) mapping done in PCI (see Figures 2, 5, and 8 in that work) and our C'^0(3 2) 
map. This filament is also suggested by the continuum e missio n maps at 1 mm and 350 fim, and 



integrated intensity HC0^(3 2) maps of 



Motte et al. 



(2003). A common correlation found in 



all these maps is that the peak integrated emission is coincident with the dust peak, where the dust 
peak is located at the refere nce position used here an d taken from our previous interferometric 



observations on this source (Cortes & Crutchei 



2006|) . However, in our N2H^ map the peak is 



offset one beam from the reference position. At this frequency, the ASTE beam is around 20" 
which coincides with the core size inferred from the interferometric data. This suggests that the 
emission from this point in the map come from independent regions of the MMl clump, which 
at the assumed distance of 5.5 kpc corresponding to 0.5 pc. Comparing with the HCO^ map, we 
see that the emission is seen enclosing a larger area around the dust peak, which includes the 
N2H'^ peak {TAv is over 45 K km s"'at the innermost contour). The integrated intensity quickly 
decreases north of the peak in both maps. In our previous maps of HCN(4 — > 3) and CS(7 — > 6) 
we saw little or almost no emission offset from the center suggesting that most of the dense 
neutral gas was unresolved and inside the MMl core. Even though CO emission is expected to be 
widespread, a similar situation is suggested in the C'^0 map where the dust peak coincides with 
a local C'^'O peak at the same position, quickly decreasing outside the dust peak. Figure [6] shows 
the N2H'^(4 — > 3) spectra towards selected offsets from the reference position corresponding to the 
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brightest points in the map. All three spectra have higher intensities than the spectra taken at the 
center of the map. Similar situations have been in seen in lo w and high mass star forming regio ns 



2007 



Zinchenko et al. 



20091 : 



Miettinen et al. 



200C 



where N^H"^ appears to be depleted towards the dust peaks (iDanie 

2, 



Johnstone et al. 



etal. 



201ol : 



2007 



Pirogoy et a 



Busquet et al. 



2011 ). It 



is likely that chemical differentiation has depleated N2H^ at the center of MMl. 



4. Discussion 

4.1. The un-identified line in the DCO+ spectrum 

To explore the nature of the un-identified line seen besides the DC0^(5 — > 4) spectrum, 
we first considered the possibility of a spurious instrumental signal. The ASTE telescope has a 
cartridge-type side-band separating the (2SB) mixer receiver working at the 350 GHz band. The 
image rejection ratio is 15 dB, which makes opposite sideband leakage unlikely. The telescope 
receiver has been carefully tested and no spurious signals are reported at that frequency (T. Sakai, 
priv. comm.). The DCO^ line was observed in the upper sideband of the ASTE receiver. Thus, 
we searched for strong lines in the lower sideband (12 GHz away), which may have leaked into 
the upper sideband. Leaking from the lower sideband would require a strong line due to the 
high side-band image rejection ratio and according to the splatalogue database for molecular 
spectroscopy, no strong lines are present in a 500 MHz interval centered at 348.169 GHz. Also, 
contamination from atmospheric lines is also ruled out as the position switching removes any 
atmospheric contribution. Another possibility is that the line is the result of a baseline ripple. 
These instabilities are usually the result of a standing wave leaked, due to reflections between 
the feed and the antenna structures, into the receiver. A possibility is oscillations of the goretex 
membrane that covers the receiver cabin under the dish. To explore this possibility, the 484 scans 
gathered were inspected one by one. Those with clear ripples were flagged-out and not used. 
After baseline removal, the remaining scans were inspected again looking for further instabilities. 
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Only those with flat baselines were kept, which imply removing about 30% of the observations. 
However even after removing spectra with ripples, it cannot be ruled out that the line feature is 
the result of such instability. Having established this, we consider the possibility of the signal 
being real. If the feature is a molecular transition, its rest frequency would be 360.1847 GHz. The 
closest molecular transition is the ^^C isotopologue of formic acid, in the 16(3, 14) — > 15(3, 13) 
at 360.18795 GHz, as reported by the splatalogue database but originally provided by CDMS 



(MuUeretal 



(lArce et al 



20051). However, the emission seen here is too strong even for '^C formic acid 



2008|) and thus, this possibility is unlikely. On the other hand, there is always the 



possibility that we are seeing a transition from an un-identified molecule. Another scenario is that 
this line corresponds to DCO^ emission at a different velocity. If that is the case, it will mean 
strong DCO^ emission from a cold component in the line of sight. This would be an exciting 



new result as deuterium chemistry is though t to be enhanced in co 



carbon bearing molecules on grain surfaces ( Bergin & Tafalla 



d cores by freezing out of 



2007|) . Nevertheless, we cannot 



conclude with certainty what the nature of the emission seen here is. It is difficult to correlate with 
emission from other molecules as their line-width is large, specially HCN. Only more independent 
observations will help confirming this un-identified emission. 



4.2. Column densities and abundances 

Assuming that the clump is in general homogeneous and isothermal, it is possible to obtain 
information about the optical depths and column densities by analyzing the line intensities and the 
velocity dispersions . Based on these assumptions, we used the radiative transfer code RADEX 



(ISchoier et al. 



20051) to estimate the column densities and abundances for several molecules 
including some of the ones presented in PCI (see Table O. The detailed source geometry is not 
known as the core is unresolved and the number of sources inside MMl is yet to be determined. 
Thus, we make no assumptions about the clump structure or possible number of sources inside. 
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As the RADEX code does not assume any kind of source geometry or velocity fields, the column 
density for the emitting material is adjusted until the line intensity is matched. We also corrected 
for beam dilution by multiplying the output from RADEX by the beam filling factor calculated 
for each transition where the source size was obtained from the interferometric observations of 



Cortes & Crutchen (|2006l) . Because the beam sizes for all molecules are comparable (differences 
are ~ 2", see Table [T]), the correction is small. Thus we should not expect strong beam-diluted 
emission if the molecular gas follows the dust distribution. However, there is evidence for 
a remarkably warm hot-core ins ide the MMl clump with T ~ 400 K and 0.03 pc in size 



( Sridharan to be submitted 



20111). As the MMl hot core is very small in size, any emission 



coming only from the hot core will be extremely diluted within any current single dish beam. 

We derived column densities for all observed molecules, with the exception of HCO^ and 

HCN as the lines are self-absorbed and more accurate geometric source modeling for the radiative 

transfer problem should be done. Instead, we used RADEX for the less abundant H^^CO^ and 

H'^CN to estimate their column densities and then used a ^^C/'^C ratio of 50 as observed in the 
II 1 

molecular ring (|Wilson & Rood, 1994 ) to derive the column densities for the main isotopomers. 
The N2H^ column density was derived using the collision rates provided by RADEX where 
hyperfine structure is considered in the calculation. However, the hyperfine structure appears 
to be blended in the main emission as the line-width is larger than the hyperfine components 
separation. The remaining molecular column densities wh ere calculated using RADEX and then 



2002bh . agreement was 



were compared with classical column density estimations dCaselli et al. 
within a factor of 2. 

To estimate abundances across the MMl main axis, we chose 3 pointings from north-east to 
south-west across MMl main axis, including the center of MMl. The pointings are indicated as 
crosses in the map (see Figure[I]). Thus, we derived column densities along the filament for 
all available molecules with the exception of DCO^ which has only 1 data point and CS which 
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does not show emission over t he 3cr level at the (30", 30 ") relative offset. The column density 



from the dust map obtained by 



Cortes & Crutched (|2006|) w as used for the cen ter pointing and 



new values were derived from the 350/zm fluxes reported by iMotte et al. 
positions. The hydrogen column densities were derived using 



(|2003|) for the remaining 



(A^H/cm-2) = 1.93 X 10^^- 



(1) 



(6 J arcsec)'^ (Z/ZQ)bT x 

where 5 y is the averaged density density flux from the source in Jy, 6^ = -\/^^^^^k^<'^^^^^ is the 
angular source size, (Z/Zq=\) is relative metalicity to solar neighborhood, is a parameter that 
reflects the variation of dust absorption cross sections (iMezger et al.lll990|) . Two values are used 
for b; b = 1.9 reproduces estimates of o-4QQfim ~ 8.3 x 10"^^ cm^(H-atom)"\ which represents 
molecular gas of moderate density ^(Hi) < 10^ cm"^, and b = 3.4 for dust around deeply 
embedded IR sources at higher densities. We used a value ofb = 3.4 in our estimate s for MMl. 



A valu e of T = 19 K is used for the dust temperature derived from the SED fitting by 



Motte et al. 



(l2003h . and x 



1.44x10** 



is the -jpj: factor for the Planck function. The hydrogen column densities 



AkT 



for the offsets are listed in Tableland also indicated as crossed in the C^^O map in Figure[Il The 
derived abundances for all molecules are listed in Table [2] and plotted in Figure |7l 



4.3. Comparison with chemical models 



4.3.1. HCN abundance 



The relative differences between molecular abundances in star forming regions can provide 
information about the evolutionary state of that region. Chemical differentiation, like the 



N/0 ratio, seems to be ubiqui tous in high mass star formation, possib 



evolutionary state of a region (IRodgers & Chamlev 



Rodgers & Chamley 



2003 



Roueff et al 



2001 



Caselli et al. 



V an indication of the 



2002bl : 



Doty et al. 



20021 : 



20071) . Besides C^^O, the most abundant molecule in 
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our sample at the center of the MMl clump appears to be HCN followed by CS, HCO"^, N2H'^, 
and DCO+. The difference be tween CS and HCN is likely within the calibration error, but larger 
for the remaining molecules. |j0rgensen et alj (120041) studied a sample of low mass and high 
mass dense cores and found the HCO"^ abundance comparable to HCN, but lower than the CS 
abundance in low-mass class I sources while in the high mass case HCN seems to be a factor of 



5 large r than HCO^ but similar to the CS abundance. A similar situation is seen by .Pirogov et al 



(120071) in a study towards a sample of southern high mass star forming regions. Due to its high 
dipole moment, the HCN molecule traces high density gas; its critical density for the 7 = 4^3 
transition is n^-nt ~ 10^ cm"^. Thus, it is expected to find HCN(4 — > 3) excited in dense clumps 
such as MMl. As the HC N molecule is most likely fo rmed in the gas phase rather than a product 



of ice mantle evaporation (IRodgers & Charnley 



20011), its abundance might indicate how evolved 



the clump is. It has been suggested that NH3 would be a major component of ice mantles and 
within the hot core region its release to the gas-phase will significantly affect the abundance of 
molecules such as HCN. As a daughter molecule, HCN is formed by dissociative recombination of 
HCNH^ which is previously formed by recombination o f C"^ with NH3. This scenario is plausible 



for higher gas temperatures (IRodgers & 



Similar models have been produced by 



Chamlev 



Doty et al 



2001) for models with T ~ 300 K in a hot core. 



20021) where they find an enhanced HCN 



abundance as the hot core evolves with time, but using H2NC'^ as the favored product to explain 
the HCN abundance. In principle the evaporation of ammonia from ice mantles will reduce 
the abundance of 0-bearing molecules and enhance the abundance of N-bearing species. Thus 
according to the previous studies, our derived abundance of HCN at the center give a possible age 
for the MMl hot core of 10"^ yr. 
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4.3.2. CS abundance 



The abundance of CS in evolved high mass clumps is likely dependent on the evaporation of 
sulphur from ice mantles by radiation from the central star. As most of the sulphur ejected from 



dust grains is li kelv locked into (iDotv et al. 



ionized sulphur (|Wakelam et al. 



20021 : 



Rodgers & Charnley 



20031), or atomic 



2004|) . this will lead to an increase of the SO2 abundance, which 



will lead to an increase in the CS abundance. Our derived value for the CS abundance towards the 
center of MMl is consistent with both models which place the CS abundance around 10"^ for an 
evolutionary time scale of lO'^ yr. However, the main issue behind using sulphur as a chemical 
clock is the uncertainty regarding the progenitors of S-bearing molecules which is still under 
debate. 



4.3.3. Deuterium fractionation 



It is known that starless low mass cores have enhanced abundances of deuterium bearing 
molecules. Due to freezing of carbon and oxygen on grain surfaces, the deute rium abundance 
can b e enhanced by orders of magnitude with respect to more evolved regions (|Bergin & Tafalla 
20071). Even though the situation in high mass clumps is not as clear, it is expected that infrared 
dark clouds will show the same increase in deuterium abundance. A useful parameter to quantify 
the importance of deuterium in clouds is the deuterium fractionation. This parameter is the 
ratio between the D-molecule and the H-molecule column densities and can be used to test the 
evolutionary stage of star forming regions. Thus, we calculated R=[DCO^]/[HCO^] = 0.0012 for 
the cent er of MMl. Fractionation values have b e en calculated for low mass class and class I 
cores by koberts & MillaJ (120001) : iTurneJ (1200 ll) : Icaselli et al. 



(I2002bl) who found values larger 



than 0.01 with most of the results about 0.06. Hi gher fractionation va 



formaldehyde, about 0.1, have been obtained by 



Bergman et al 



In the case of high mass cores, results are not as numerous. 



ues based on deuterated 



(1201 11) towar d .s p-O phiuchi 



Miettinen et al. 



(20091) found 
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moderate deuteration of about 0.03 to 0.04 for Orion B9; while 



Chen et al 



(1201 oh found moderate 



deuteration between 0.01 to 0.05 toward the infrared dark cloud G28. 34+0.06. 



Roueff et al 



(|2007|) modeled deuterium fractionation in warm dense clumps, such as MMl, and found that 
[DCO"*"]/[HCO"^] decreases rapidly with temperature giving fractionation values similar to our 
findings for T < 30 K, which suggests that MMl might be more evolved than previously thought. 
Considering that infall is proceeding, as showed by PCI, and a hot core has already formed, it is 
possible that gas in the MMl envelope has been heated by the hot core to the point where the 
dust and the gas are no longer in LTE. At the same time, it is also possible that the dust emission 
is optically thick. In both cases the dust temperature is not representative of the real kinetic gas 
temperature inside the MMl clump. Even though the amount of deuteration found in MMl is not 
significant as in pre-stellar sources, it is still significant for a warm environment such as the one 
suggested by the model. The DCO+ molecule can be assembled by combinati on of CO with H^D "^ 
which is very efficient at low temperatures. However, chemical modeling by Roberts & Millar 
(|2000|) suggests that this mechanism is n o longer efficient when temperature increases, which 



is also suggested by 



Roueff et al. 



(|2007|) . As the envelope gas temperature rises the relative 
abundance of C, N, and O molecules tend to decrease the abundance of H2D+ which can no longer 
efficiently combine with CO to produce DCO^. Instead, CH2D+ might become the dominant 
pathway for DCO^ creation by CH2D+ + O — > DCO^ + H2 as the deuteration of CH3 can proceed 
to much higher levels than H2D'^, the usual way to produce DCO^ (A. Wootten, priv. comm.). As 
previously mentioned, we question how well the derived temperature from the SED represents the 
real envelope gas temperature; it could well be warmer than the dust as the dust may have become 
optically thick. So it is possible that the DCO^ emission that we are seeing comes from a different 
chemical pathway than what is seen in colder sources. 
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4.3.4. N2II* and HCO* abundances 



For high mass star forming regions, HCO"^ abundances hav e been derived from a 



sources similar to MMl with values in the 10 ''to 10 '° range (iZinchenko et al. 



sample of 



2009). 



column densities have been derived as well, toward several h igh mass cores at different stages o f 



evolu tion with values that range between lO''^ to 10'^ cm (|Szymczak et al. 



2007 



Sakai et al. 



2010h . Our derived HCO^ abundance towards the center of the MMl clump is X(HCO^) = 
1.5 X 10"'' and the H'^'CO^ column density is 2 x 10'^ which are consistent with results in 
similar regions and with the chemical models previously described. The N2H^ abundance has 
also been derived towards high mass cores, mostly from N2H^(1 0) observations as there are 
few N 9H^(4 — » 3) results in the literature. Current findings vield abundances between 10~ " to 



10"'" (Daniel et al. 



20071 : 



Pirogov et al 



20071 : 



Zinchenko et al 



2009 



Miettinen et al. 



2009). Our 



derived abundan ce towards the center of MMl is X(N2H"^) = 9.6 x 10 " which is consistent with 

(|2011r ) produced a chemical model for the high mass star forming 



previous results. iBusquet et al. 



region AFGL 5142. They modeled the effect of a central source on the abundance of several 
species including N2H^. At an evolutionary time of 10'' yr, their model predicts an abundance of 
X(N2H"^) ~ 10"'° for the central core of AFGL 5142, which is about the abundance value that we 
find for MMl and in the same order of magnitude with the age obtained previously by comparing 



Lintott et al 



(|2005|) compared 



HCN, CS, and DCO^ abundances with other chemical models, 
the production of N2H'^ and related chemical species for high mass star formation environments. 
While they do not present abundances for N2H^, they studied the relative abundance between CS 
and N2H^ which we can compare here. Their model predicts that accelerated infall will produce 
a differentiation between CS, HCN, and N2H^, where they expect that relative abundances will 
decrease due to the destruction of N2H'^ by CO and enhancement of CS by the ac celerated 



Lintott et al. 



collap se. At the center of MMl we obtained a log([CS]/[N2H+]) = 2.3, which in the 
(|2005|) model is close to no acceleration scenario. As MMl is infalling with speed slightly 
supersonic but sub-alfvenic (see section 6.6), it is difficult to reconcile this with a no acceleration 
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model. However, we do see an enhancement of the CS abundance at the center of MMl which 
will be discussed in the following section. 
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Table 2: Column densities and abundances are presented here. 



Molecule 


Offsets 


Column Density 


Abundance 


Line- width 




[(arcsec, arcsec)] 


[cm-^] 




kms 


DCO+ 


(0,0) 


1.2 X 10^^ 


8.8 X 10-13 


2.9 


N2H+ 


(0,0) 


1.3 X lO''* 


9.6x 10 " 


7.3 


N2H+ 


(30,30) 


6.0 X 1013 


1.1 X 10-1^ 


6.3 


N2H+ 


(-30,-15) 


1.3 X lO^"^ 


2.0 X IQ-^^ 


5.6 


HCO+ 


(0,0) 


1.0 X 10^^ 


7.4 X 10-^^ 


8.4 


HCO+ 


(30,30) 


2.5 X 10'^ 


7.9 X 10""' 


6.2 


HCO+ 


(-30,-15) 


1.5 X 10^^ 


2.4 X 10-^^ 


6.2 




(0,0) 


2.0 X 10^3 


1.5 X 10-11 


3.0 


HCN 


(0,0) 


2.0 X 10^^ 


1.5 X 10"^ 


9.5 


HCN 


(30,30) 


8.0 X 10^^ 


2.5 X 10-'' 


7.6 


HCN 


(-30,-15) 


2.5 X 10^^ 


4.0 X 10-^ 


14 


Ri^CN 


(0,0) 


4.0 X 10^^ 


2.9 X lO-i'' 


4.8 


CS 


(0,0) 


1.0 X 10'^ 


7.3 X lO-'' 


4.3 


CS 


(-30,-15) 


1.3 xlO^^ 


2.1 X lO-'' 


9.4 


C'^O 


(0,0) 


2.9 X 10^^ 


2.1 X 10-^ 


6.4 




(30,30) 


8.0 X 10^5 


2.5 X 10-^ 


5.1 




(-30,-15) 


1.1 X 10^^ 


1.7 X 10-^ 


5.7 
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Table 3: Hydrogen column densities for three selected pointings are shown here. 

Offsets (or, 6) Column Density 
([arcsec], [arcsec]) [cm~^] 
(-30,-15) 6.3 X 10^3 

(0,0) 1.5 X 10^4 

(30,30) 3.2 X 10^3 
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4.4. Spatial differentiation 



In all the maps prod uced from our obse rvations, we see evidence for a filament in the MMl 



clump (see also PCI and 



Motte et al. 



(I2003|) ). The C^^O integrated emission map is consistent 
with this, as well as our abundance derivation which is constant along the main axis of the 
filament. The spatial distribution of both '^CO and dust (both 850 and 350 /im) also show a 
single filament with a continuous progression in the emission from the MMl center to the end of 
the filament at the south-east edge. The '^CO map however, shows strong peaks in MMl and at 
the south-east edge where there is a likely boundary with the H II region. Only the MMl center 
peak is seen in dust emission maps. Our C^^O map shows two distinct areas with strong emission 
around the dust peak and close to the south-east boundary. There is a clear region in the middle 
where the emission drops in intensity, which is not seen in the dust maps and it is not completely 
clear in the '^CO map. This suggests that the MMl filament might not extend completely to the 
boundary with the H II region as C'*^0 is likely showing us the integrated emission through the 
cloud rather than just the envelope. The influence of H II region ionization and shock fronts as 
external pressure over the MMl clumps is important as a mechanism for triggered star formation. 
The N2H^ emission is also distributed like C^^O and HCO^, but we noticed that the peak emission 
in N2H"^ is not coincident with the dust peak as seen with the other molecules. The main peak in 
the N2H^ integrated emission is 1 beam away from the ce nter to (0", 20"). Thu s, it cannot be 



Caselli et al 



(|2002b() suggested that 



attributed to an oversampled grid re-distributing emission. 
N2H^ followed the dust emission in cold cores as it has been seen that N-bearing species do not 
get fr ozen on dust grains unt il densities reach the 10^ to 10^ range with gas temperatures around 



10 K (IBergin & Tafalla 



20071). However, when the gas gets warmer due to heating from the central 



source(s), CO is rapidly released from grain ice mantles quickly destroying N2H^. If the gas in 
MMl is being warmed by the hot core releasing CO into the gas phase, it may explain the offset 
seen in the N2H^ peak emission with respect to the dust peak. As our sub-millimeter observations 
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trace higher critical densities, it is likely that we are seeing gas closer to the hot core where the 
temperatures are likely higher than the dust SED suggested. 

We also noticed that the emission from the neutral molecules is more compact than the 
ions, with the exception of the H'^CO^ and H^^CN which only have significant emission from 
the center of MMl. The HCN and CS emission are mostly confined to the clump, where CS is 
remarkably compact (see PCI Figures 7 and 8). But, the N2H^ and HCO^ integrated emission 
are much more widespread over the filament. All these molecules have similar critical densities 
for the rotational transitions observed; thus, they presumably come from gas under the same 
physical conditions. In order to explore possible explanations for this spatial differentiation, we 
will consider the effect that the magnetic field threading MMl might have over the gas along the 
filament. If the field is frozen in the ions, it is likely that the tension created by the field will limit 
the movement of ions along the filament; while the neutrals can diffuse more freely through the 
ions. Table [2] shows the line- widths obtained by adjusting Gaussian profiles to the line emission 
at the selected positions along the strip. It is evident that the line- width of N2H'^ and HCO^ are 
consistently smaller than the neutral molecule line- widths. Particularly at the center we also have 
DCO^ and the ^^C isotopomers of HCN and HCO^. As these species are optically thin their 
velocity dispersions are more representative than their more abundant ^^C isotopomers. All ions 
show smaller velocity dispersions than neutrals with DCO^ and H^^CO^ about 3.0 kms~'FWHM ; 



while H CN has the largest line-width at 4.8 kms . It has already been noted by 



(12000, 



I 



Houde et al 



lines. 



[b|) that ions s hould have smaller velocity dispersion due to trapping on magnetic field 



Li et al. 



(|2010h proposed that the smaller v elocity dispers ion of ions with respect to neutrals 



is a signature for turbulent ambipolar diffusion (|Li et al 



20101). While ambipolar diffusion is 



inevitable at som e point during the clump evolution, how fast the neutrals will diffuse is still a 



matter of debate (McKee & Ostrikei 



2007 ) 



Chemical models of dense clumps, including inner hot cores, have calculated radial 
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profiles for the abundan ce of different molecules (|Doty et al. 



Nomura & Millar 



2OO2I : 



Rodgers & Chamley 



2003 



2004|) . All these models find that ion abundances increase rapidly with distance 
from the central source(s). Unfortunately, not all these models present results for all the molecules 
presented here. Thus, we will con centrate on comparing with the model that presents the largest 



number of molecular abundances. 



Rodgers & Chamleyl (|2003[ ) produced models for envelopes 



surrounding hot cores, as we believe is the MMl case. They developed models based on the 
dynamical state of the envelope, considering static and collapsing cases (unfortunately without 
including magnetic fields). For their collapse model and assuming gas temperatures of 300-400 K 
for the hot core, they found HCO^ abundances up to a few tenths of 10"^ at 1 pc from the central 
source(s) for 10"^ yr and 10"^ for 10^ yr of evolution after collapse started. As PCI presented 
compelling evidence of infall in MMl, comparing with this model seems pertinent. Even though 
the UMIST database includes NiH"^, they did not present predictions for that molecule in their 
work. Table [2] shows HCO^ column densities and abundances at distances of about 1 pc from 
the center at both sides of the filament. Our values are in the range of 10 ~'° to 10"' ' which is 



Rodgers & Charnleyl (l2003h model. The HCO 



at least an order of magnitude lower than the 
abundance is derived from H'^CO^, which is moderately optically thin (from the RADEX results, 
T = 0.7), but we assumed a '^CO/'^CO ratio of 50 which is not necessarily accurate. However 
even using values close to 70 fo r the '^CO/'^CO ratio, the H CO^ abundance will not agree with 



the prediction. Interestingly the 



Rodgers & Chamleyl (|2003h model predicts a sudden decrease in 



the abundance of neutral molecules such as CS, with a cutoff point around the collapse front. The 
case of HCN is not as clear, as their models show a small decrease in abundance at the collapse 
front, but an increase towards the envelope. The HCN abundance predicted at 1 pc from the center 
is about 5 x 10"^, which is about what we see from our data. For CS we derived an abundance 
similar to HCN at 1 pc, but the situation is not symmetric as we do not have emission at the other 
end. This CS emission at offset (-30", -15") is likely coming from the same place where we see 
additional HCN emission ( see Figure 7 in PCI); thus, it would not be a valid comparison with 
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the model. However, we do know that th e CS emission in MMl is fairly compact which is also 



predicted by 



Rodgers & Chamley 



(l2003h . The overall picture presented by the model is similar 



to what we have determined from our observations, particularly with the strong CS cutoff at the 
collapse front. Additionally, the panel presented by PCI for HCN infalling spectra places the 



Rodgers & Chamley 



collap se front at 25" in radius (or about 0.9 pc ) which is also similar to the 
(|2003h results. Even tough our abundances for neutral molecules agree with that model, the HCO+ 
molecule does not. However, it is feasible that changes in the chemical abundances due to the 



heating of the hot core explains the spatial differentiation that we are seeing in the MMl clump. 



4.5. The ionization rate 



4.5.1. The electron fraction 



The abundance of molecular ions depends on a complicated interplay between gas-phase 
chemistry and cosmic -ray driven chemistry. Molecular ion abundances are critical to understand 
the interaction between the magnetic field and the neutral gas as the magnetic field can only 
influence the dynamics of dense neutral gas through the ions. The electron fraction is estimated 
by charge equilibrium between the electrons and the most abundant ions in a molecular clump 
as we assume overall charge neutrality. The electron fraction permits calculations of neutral-ion 
coupling parameters, which we can use to understand the impact of the magnetic field in the 
gas dynamics of the infalling MMl clump. As H+, H7D^, DCO^, HCQ-", N7H+ , and N7D+ ar e 



considered the most abundant ions in a high mass clump (iMiettinen et al. 



2008L 



2002b : 



Bergin et al. 



19991 : 



Williams et al. 



1998 



Wootten et al 



20091 : 



Caselli et al. 



1979h . a lower limit for X(e) 



can be written as 



X(e) > X(H+) + X(H2D+) + X(DCO+) + X(HCO+) + X(N2H+) + X(N2D+) 



(2) 
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As previously shown, the deuterium fractionation is not large in MMl as the clump has already 
evolved from a pre-stellar phase. This suggests that deuterium species, such as N2D"^ and 
DCO^, are not directly important in the derivation of the electron fraction and can be neglected 
in the calculation. The case of HiD^ is somewhat different, as it is directly related to the H3 
abundance through H3 + HD <r-> HiD^ + H2, and thus ca nnot be negle c ted fr om equation [21 The 



abundance of H2D^ can be estimated from the survey of Caselli et al. 



(l2008h towards a sample 



of star forming dense cores. In that work, the abundance of H2D"^ was calculated by observing 
the ortho-H2D+(li,o - li,i) line with CSoE. To obtain an estimate of the H2D^ abundance, we 
used an average of their values towar d dense cores, or X(H7 D^) = 2.3 x 10~"', and assumed an 



ortho-to-para ratio of 10 as used by 



of the abundance values derived by 



Miettinen et al 



Caselli et al 



j2009l) . We are here assuming an average 



(|2008b are representative for MMl. Some 



through the relation derived by 
simplified c 



dense cores used by them are indeed similar to MMl (such as Ori B9), so using an average 
is a compromise for a better representation of H2D^ abundance in star forming regions. Now, 
the deuterium fraction R provi des a rough estimat e of the r = [H2D^]/[H3 ] abundance ratio 

ipsi et all (120041 ). R = {r + 2r^)/(3 + 2r + r^), by assuming a 
l emica l network for deuterium in molecular clouds. Thus, by numerically solving the 
(120041) relation and using = [DCO+]/[HCO+] = 1.2 x lO""* from our observations, 
we obtain ^(Hp = 6.4 x 10 *^ giving an electron fraction of X(e) = 6.5 x 10"^. Our estimation 
for the electron fractio n is consistent with other high mass star fo rming sites such as O ri 9 where 



Crapsi et al. 



1 r 



Miettinen et al., (,20091) calculated 6x10"^ and DR21(0H) where 



Hezareh et al 



3.2 X 10 ^, and in the upper range of low mass cores such as L1544 where 
found X(e) ~ lO^^ 



(l2008b obtained 



Caselli et al. 



(l2002bl) 



■*The CSO has a similar beam size as our ASTE observations 
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4.6. The Ion-Neutral coupling 



The effect of a magnetic field in a weakly ionized gas depends on how well the neutrals are 
coupled to the ions. A well coupled field-gas will be more efficient in restraining the neutrals 
movement in directions parallel to the field where magnetic support does not directly delay the 
gas infall. We first noticed that infall in MMl is slightly supersonic but sub-alfvenic, where 
Vif = 0.5 kms^Hs the infalling speed, Cs = 0.46 kms"Hs the sound speed, and Va = 2.4 kms"Hs 
the Alfven speed. These parameters yield a Mach number of M=l.l and Alfvenic Mach number 
of Ma = 0.2. As the infall is sub-alfvenic, MHD waves can in principle propagate faster than the 
collapsing gas. If the gas is warmer than the dust as previously suggested, the sound speed will be 
larger making the infall sub-sonic. 

It has been found in laboratory plasmas that magnetic Reynolds numbers much larger than 
1 suggests that coupling between the gas and the field is strong and in principle, this allows the 
gas and the field to move as a single system. This number is defined as the ratio between the 
inertial term to the viscous term in a magnetized fluid, or R,„ = where Vif is the fluid speed, / 
is the length scale in which the magnetic field has a significant change in magnitude, and vm is 
the magnetic viscosity defined as vm = Va^/v„,, where v,„ is the ion-neutral collision frequency 
defined as y,„ = X(e)nH2 < cr,„v >, where < cr,„v > is the ion-neutral collision rate. Thus, we 
calculate the magnetic (ambipolar diffusion) Reynolds number as Rm = ^X(e)nH < o-jnV >, 
w here / is calculated t o be a bout ~ 5" or ~ 0.1 pc obtained from the dust polarization map 



of 



Cortes & Cmtcheil (120061) . The Alfven speed is Va = 2.2 x l O^B/(n„A„Y'^ where we take 



An = 1 .4 f or representative mo 



cm from 



ecular gas of cosmic abundance (|Zweibel 



20021), tiH = 4.5 X 10^ 



Cortes & Crutcheii (l2006h . and B = 855 fiG from PCI. The < cr,„v > is the ion-neutral 



collision rate which is calc ulated using the Langeving approximation for the polarization potential 



(IMcDaniel & Mason 



19731). Here, we use a value of < cr,„v > 



gas in a typical star forming region (IMouschovias & Paleologou 



.69 X 10 ^ cm"* s ^ for molecular 



19811) • Thus by using our derived 
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values for X(e), we calculate = 18.1. As mentioned before, the ambipolar diffusion Reynolds 
number describes the stability of a magnetized flow, or how close the fluid is to become turbulent. 
The boundary, or critical point, is often found for astr ophysical fluid plasmas comparable to 



the Hartmann number M ( Myers & 



■Giersonskv 
idsoJ 



19951) . which is defined as the ratio of the 



200 lb and it can be written sls M = Bl sJcr/T], 



Lorentz force to the viscous forces (iDavidsor 
where cr = c^KAttvm) is the electrical conductivity a nd rj = pv is the dynamic viscosity function 



of the kinematic viscosity v. 



Myers & KhersonskyI (|1995|) derived the Hartmann number as a 



function of column density as M = 1.99 x 10^(A'^/10^^cm"^)X(e)'^^, the expression that we use 
here. In astrophysical weakly ionized plasmas, if 7?^ » M » 1 the flow and the field can be 
considered to be turbulent as the dissipative term is not large enough to damp the inertial term, 
but if M > i?m » 1 the field can still damp the turbulence in the flow stopping the cascade into 
smaller eddies. Thus, we directly calculate the Hartmann nurn ber finding M = 7 x 10^. This value 



is consistent with calculations of 



Myers & Khersonskvi (i l995i) towards similar dense cores and to 



recent numerical simulations of sub-alfvenic flows ( Li et al 



2008 



McKee et al. 



201 



The ratio between the core size, R, and the minimum MHD wavelength, or cutoff wavelenght, 
for propagation of MHD waves has also been used to characterize the ion-neutral coupling. This 
ratio is known as the wave coupling number W = R/Ag, where '^a = Va is the Alfven speed . 



and Vi n is the ion -neutra' 



1995 



Bergin et al 



collision frequency (lElmegreen & Fiebig 



1993 



Myers & Khersonsky 



1999|) . Since free MHD waves propagation in a clump requires minimum 
wave damping, we expect this minimum wavelength to be less than the observed length-scale, 
orW> 1, as this ratio giv es the range between the maximum possible \yavelen gth (core size) 



and the cutoff wavelength (Kulsrud & Pearce 



19691 : 



Elmegreen & Fiebig 



1993b. Thus, when 



W » 1 we expect that the field and the gas are coupled over a large range of length- scales. 
Using our previous determination for the strength of the magnetic field in this clump, the Alfven 
speed, and the ion-neutral collision frequency we obtain a value for the wave coupling number 
ofW= 110. Previous determinations of W suggest that values close 1 are obtained when the 
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field is in equipartition with self-gravity; while greater values are consistent with a dominant 
magnetic field and strong coupling bet ween the field and the neutr al gas. Values around W ~ 10^ 



have been found toward dense cores (IMyers & Khersonskv 



995). but smaller values have 



1999b. Values smaller than 10^ 



also been obtained (W ~ 20) for massive cores (iBergin et al.l 
suggest marginal coupling, where MHD waves can still propagate, but damping by the field is 
still possible. We have obtained a result larger than the equipartition value, but still small in 
comparison with other dense cores. As we have previously found, the mass-to-magnetic-flux ratio 
in this clump is super c ritical and the dispersion f unction in the field lines, derived by PCI from 



the polarization map of 



Cortes & Crutchen ilOOff ). suggests that a significant fraction of the field 



is turbulent. Thus, a sligthly larger than 1 value for the ambipolar difi'usion Reynolds number and 
a larger than 1 wave coupling number value suggest that in MMl the neutral gas has already or is 
in the process of decoupling from the field and a transition between a more dominant field and the 
full onset of gravitational collapse is ongoing. However, there is still sufficient dynamic range in 
the length-scales for the gas and the field to move as a single system. It is interesting that we are 
seeing an infalling clump with a field just decoupling from the gas and a hot core already formed 
where molecular abundances are being changed by the central source(s) and where no powerful 
outflows have been confirmed yet. This points to an intermediate scenario for the evolution of 
a high mass star forming region. The clump seems to have evolved creating the inner source(s) 
before the envelope stops accreting and before the gas gets fully decoupled from the ambient 
magnetic field. Moreover and even though MMl is collapsing, we see evidence for an envelope in 
which the ions have a more extended distribution than the neutrals from species excited at similar 
critical densities. In principle ambipolar difussion can distribute the magnetic flux to the envelope 
restricting the movement of ions while the neutrals collapse. This mechanism can operate in the 
presence of MHD turbulence in which numerical simulations suggests that the ambipola r diffusio n 



time scale is significantly reduced, but it can also efficiently redistribute magnetic flux (|Li et al. 



2008h . Nevertheless, we cannot rule out that the spatial differentiation seen here is produced by 
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chemistry alone, we favor the previous mechanism as the more important one where chemical 
differentiation is primarily responsible for the spatial differentiation seen in MMl. 



5. Summary and Conclusions 



We have presented new and re-processed sub-millimeter data from W43-MM1. New 
N2H+(4 3), DC0'"(5 4), and SiO(8 7) data have been obtained as well as previous 
un-published C'*^0(3 —> 2) observations. We have re-processed previous HC0^(4 — > 3), 
Hi3co+(4 ^ 3), HCN(4 ^ 3), H'3cn(4 ^ 3), and CS(7 ^ 6) data to obtain chemical 
abundances in MMl. New un-published re sults with the SMA, place a war m hot core, of about 



400 K, inside MMl as well as an outflow (ISridharan to be submitted 



20111). We did not detect 



SiO emission over the 3cr threshold; thus, the outflow emission seen in the SMA data is not 
seen from single dish measurements as it is likely deep inside the clump or its emission is 
highly diluted within the single dish beam. We have calculated column densities by using the 
RADEX code, which agree with classical methods within a factor of 2. Abundances for these 
lines have been derived for 3 selected points alo ng the main axis of the MMl clump. W e have 



compared our findings w i th chemical models 



Nomura & Millail (120041) 



RouefF et al 



by 



Doty et al 



(120021) : 



Rodgers & CharnlevI (l2003h 



(|2007h . Comparison with these models put the MMl 
clump hot core within an evolutionary age of about 10"^ yr. We also found the N2H"^(4 — > 3) peak 
emission offset from the dust peak by about 0.5 pc, which can be explained by CO enhancement 
at the MMl center due to the hot core as CO quickly destroys N2H^. We derived a deuterium 
fractionation of 0.0012 based on our [DCO^]/[HCO"^] abundance ratio. This fractionation is 
consistent with an already evolved clump in which the hot core is releasing C and O bearing 
molecules altering the chemical composition of the MMl clump and is also consistent with the 
depletion of N2H^ seen at the dust peak. These findings question how well the dust temperature 



is coupled to the gas temperature, as comparison with models suggest higher gas temperatures 
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than the 19 K derived from the dust SED. We conclude that the gas should be warmer than the 
dust inside the clump as the dust emission might be likely optically thick. At the same time, our 
maps of integrated emission suggest different spatial distributions for ions and neutral species; 
while the ions are more wide-spread around the MMl center, the neutrals (except C^*^0) show a 
compact ditribution around the dust peak. It is possible that magnetic tension is restricting the ion 
movement while the neutrals collapse more freely, but we favor chemical differentation induced 
by the hot core as the main mechanism behind the spatial differentation seen in MMl. The 
map suggests that the MMl clump does not follow completely the filament suggested by the dust 
and the ^^CO maps. Two regions are seen, one corresponds to the MMl clump and the other to 
what we believe is the interface with the H II region (a PDR?). If this is confirmed, then the shock 
front from the H II region has not yet reached the MMl clump, which may have implications for 
triggered star formation models. 

By using our derived abundances for the ions, we have derived an electron fraction of 



X(e) = 6.5 X 10 ^ where we used an average value for the H2D^ abundance from 



Caselli et al. 



(|2008r) and our deuterium fraction to derive the H3 abundance. Our value for X{e) is consistent 
with similar values in other high mass star forming regions found in the literature. Using the PCI 
estimation for the magnetic field strength, we calculated an ambipolar diffusion Reynolds number 

= 18.1, a Hartmann number M = 7 x 10^, and a wave coupling number W = 110. These 
values suggest that the field has or is in the process of decoupling from the fluid as infall proceeds. 
As the Reynolds number is still not comparable with the Hartmann number, the field may still 
damp MHD waves, but as suggested by the wave coupling number, there is still enough dynamic 
range in the length-scales for the field and the fluid to move as a single system. By combining 
our findings here, we believe that the MMl clump is in an intermediate stage of evolution. At 
an intermediate stage of evolution, the inner source(s) has(have) formed a hot core where a 
strong outflow is detected. The envelope is still contracting and it is not fully decoupled from 
the ambient magnetic field. To explore higher resolution spatial scales, particularly at the cutoff 
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(minimum wavelength) level, interferometric observations will be needed. Thus, ALMA is the 
perfect instrument to study W43 further. 

The National Radio Astronomy Observatory is a facility of the National Science Foundation 
operated under cooperative agreement by Associated Universities, Inc. 
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Fig. 1. — A 60" X 60" integrated intensity map of C'^0(3 2) from W43-MM1 centered at 
the dust peak. The color scale correspond to the integrated emission after applying a nearest 
interpolation scheme. The contours are choosen to be 8, 12, 16, 20, and 24 K km s^^ while the 
colorscale is indicated in the colorbar. Crosses indicate the positions at which we derived fractional 
abundances from our sample of molecules. The offsets are (30", 30"), (0", 0"), and (-30", 15") in 
{a, 6). Note tha the reference position is offset from the center of the map by 60" in a and 30" in 
6. 
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Fig. 2. — The figure shows our SiO(8 — > 7) spectrum towards W43-MM1. The black line marks 
the zero baseline and the systemic velocity of the clump, 98.9 kms Hs shown by the green line. 
The spectra is binned every 10 channels giving a resolution in velocity of about 1 km s^'giving and 
rms noise of 0.24 K. The intensity upper limit at the systemic velocity is 0.6 K. 
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Fig. 3. — The spectrum shows our DC0"^(5 — > 4) detection towards W43-MM1. The Gaussian 
fit is shown by the red curve, note that the systemic velocity of the clump, 98.9 km s^^is shown by 
the green line. Also, the un-identified line is also shown with its corresponding Gaussian fit (blue 
curve). 
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Fig. 4. — A 60" X 60" integrated intensity map of N2H"^(4 3) from W43-MM1 centered at the 
dust peak is shown. The color scale correspond to the integrated emission after applying a nearest 
interpolation scheme in units of K km s"^. The contours are chosen to be 2.5, 5, 7.5, 10, 12.5, and 
15 Kkms ^ 
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Fig. 5. — The spectrum shows the N2H"^(4 —> 3) from the center of W43-MM1. The Gaussian fit 
is shown by the red curve, note that the systemic velocity of the clump, 98.9 kms Hs shown by 
the green line. The spectra is binned every 10 channels giving a resolution in velocity of about 1 
km s^^ The spectrum was obtained with a higher signal to noise than the data presented in Figure 
IHas we integrated additional time at the center of MMl. 
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Fig. 6. — Spectra from N2H^(4 — > 3) towards selected pointings in W43-MM1. The offsets are 
respect to the main reference position an indicated by the label box in each spectrum. The spectra 
is binned every 10 channels giving a resolution in velocity of about 1 km s"' . 
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Fig. 7. — The plot shows the fractional abundances derived for all the species observed at the 
three selected pointings along the MMl strip (main axis). The points are indicated by their relative 
distance respect to the center. Thus, (30", 30") is at 1.1 pc, while (-30, -15) is at -0.9 pc from the 
center. Each molecule is represented by a symbol indicated in upper left box. 
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